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Abstract 

We consider a stabilized RSI model in the energy range below the direct produc- 
tion of KK states. In this range we work out the effective Lagrangian due to exchange 
of heavy KK tensor graviton and scalar radion states and compute explicitly the cor- 
responding effective coupling constants. As an example, the Drell-Yan lepton pair 
production at the Tevatron and the LHC is analyzed in two situations, when the first 
KK resonance is too heavy to be directly detected at the colliders, and when the first 
KK resonance is visible but other states are still too heavy. It is shown that in both 
cases the contribution from the KK invisible tower leads to a modification of final par- 
ticles distributions. In particular, for the second case a nontrivial interference between 
the first KK mode and the rest KK tower takes place. Expected 95 % CL limits for 
model parameters for the Tevatron and the LHC are given. In the Appendix useful 
formulas for the cross sections and distributions of various new 2 — > 2 processes via 
heavy KK tower exchange are presented, the new formulas containing nonzero parti- 
cle masses for final state fermions and bosons. The formulas and numerical results 
are obtained by means of the CompHEP code, in which all new effective interactions 
are implemented providing a tool for simulation of corresponding events and a more 
detailed analysis. 

1 Introduction 

The paradigm of modern quantum field theory implies that the fundamental short-range 
interactions of elementary particles are either mediated by massive bosons or involve con- 
finement. In particular, the weak interactions are mediated by massive gauge bosons. 

It is a common knowledge that historically the weak interactions were first described 
by contact four-fermion interactions, because the large masses of the intermediate gauge 
bosons disguised the nonlocal character of this interaction for small energy or momentum 
transfer. Nevertheless, this approach enabled the physicists to treat the weak interactions 
theoretically long before the consistent theory was formulated and the intermediate gauge 
bosons were discovered. 

Nowadays, as the high energy physics is looking for new interactions beyond the Standard 
Model (SM), it is very likely that we are in a situation, which is reminiscent of Fermi's theory 
of the weak interaction. 

There are many theoretical schemes, which predict new interactions beyond the Standard 
Model, mediated by new particles, but these new particles may be too heavy to be directly 
found in experiments. Thus, it is worthwhile to consider the situation, where the energies ac- 
cessible at the existing and the upcoming colliders are well below the threshold of production 
of these new particles. In this case the new interactions, predicted by a particular model. 



are reduced to contact interactions of the Standard Model particles, which are defined by 
the model at hand. 

The usage of contact interactions, or effective higher dimensional operators, is a well- 
known way of parameterizing possible deviations from the SM in a model independent way 
[H |2l [3l H] . Such operators are introduced by hand in phenomenological extensions of the 
Standard Model, the only restriction on their form usually being the conservation of the 
Standard Model symmetries. In various processes at hadron and lepton colliders, the effective 
operators may be probed with the aim to get some first indications of a manifestation of 
physics beyond the SM or to obtain restrictions on the parameters of the effective Lagrangians 
(see, e.g., [SI [6l [TJ [8] ) . However, there appears a large number of admissible effective operators 
each coming with its own coupling constant, which results in an uneasy problem of extracting 
many parameters from the experimental data. 

Consideration of models with extra dimensions (see, [U dUl [HI [121 [131 [HI [H]) leads 
to a very definite prediction for the structure of the contact interaction operators entering 
the effective Lagrangian. In particular, the contact interactions arising in such models are 
universal in the sense that they are characterized by only one dimensional constant. An 
experimental observation of such contact interactions could be a strong argument in favor 
of models with extra dimensions. Contact interactions due to summation of the exchange 
of Kaluza-Klein (KK) towers within the ADD (Arkani-Hamed, Dimopoulos, Dvali) scenario 
were studied in [16]. The collider phenomenology of the contact interactions appearing 
below the production threshold of KK modes in the RSI model, such as changes of dis- 
tribution tails, was discussed in [17] . Contact interactions were also considered in theories 
with warped universal extra dimensions, where contributions of KK vector boson towers to 
Fermi's constant were estimated [13] . 

In this paper we study the stabilized RSI (Randall, Sundrum 1) model below the pro- 
duction threshold and perform a more accurate derivation of the effective contact interaction 
Lagrangian. The latter enables us to take into account the interactions of the scalar compo- 
nent of the multidimensional gravity and to demonstrate explicitly that its contributions to 
the contact interaction are much smaller than those of the tensor modes. From the analysis 
of Drell-Yan distribution tails including KK contributions and those of the SM with all the 
modern uncertainties and a natural restriction on the KK resonance width to be smaller than 
its mass, we give expected collider bounds on the effective Lagrangian parameters for both 
the Tevatron and the LHC. For the case where the mass of the first KK tensor mode lies 
within the collider energy reach we show how the Breit-Wigner distribution of the resonance 
is modified due to the contributions of all the remaining modes including the destructive 
interference with the resonance. 

The paper is organized as follows. First, within a stabilized Randall-Sundrum model, we 
derive an effective Lagrangian for the interactions of the Standard Model particles induced by 
Kaluza-Klein excitations in the case, where the center of mass energy is below the threshold 
of the excitation production. An important point is that we explicitly calculate the effective 
coupling constants. Next we discuss some collider manifestations of the effective contact 
interactions giving the relevant formulas in the Appendices. In particular, the formulas for 
the cross sections and distributions include both tensor and scalar contributions and take 
into account the masses of the final state particles. 
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2 The Effective Lagrangian for the stabilized RSI 
model 



The characteristic feature of theories with compact extra dimensions is the presence of towers 
of Kaluza-Klein excitations of the bulk fields, all the excitations of a bulk field having the 
same type of coupling to the fields of the Standard Model. If we consider such a theory for the 
energy or momentum transfer much smaller, than the masses of the KK excitations, we can 
pass to the effective "low-energy" theory, which can be obtained by the standard procedure. 
Namely, we have to drop the momentum dependence in the propagators of the heavy modes 
and integrate them out in the functional integral built with the action of the theory. This 
can be easily done, if the self-interaction of the modes is weak, and one can drop it as well. 
As a result, we get a certain contact interaction of the Standard Model fields for each bulk 
field of the multidimensional theory. If we also assume that the fundamental energy scale 
of the (4 + (i)-dimensional theory M is of the same order of magnitude, as the inverse size 
of extra dimensions, then the masses of the KK excitations are proportional to this energy 
scale M, and the wave functions of the modes are proportional to M^^'^. This defines the 
coupling constant of the contact interactions up to a dimensionless factor. The particular 
structure of the contact interaction Lagrangian is fixed by the corresponding structure of the 
SM current coupled to the zero mode of a bulk field and by the spin-density matrix of its 
KK modes. This leads to a number of very concrete predictions for collider phenomenology. 

The bulk field, which appears in any theory with extra dimensions, is the gravitational 
field. The question, what fields, besides the gravitational one, propagate in extra dimensions, 
has no unique answer yet. The theory of universal extra dimensions [T9] allows all the fields 
of the Standard Model to propagate in the bulk. Other approaches allow only some of the 
Standard Model fields to propagate in extra dimensions. In particular, another interesting 
assumption is that only the gauge fields live in extra dimensions. A motivation for this can 
be the theory of "fat" branes, where the fields of the Standard Model are trapped on the 
brane by a particular mechanism. It turns out that it is easy to trap the fermion fields on 
the brane, but there is no mechanism yet for trapping the gauge fields |20j . 

There are two main approaches in theories with large extra dimensions, - the ADD 
scenario and the Randall- Sundrum model, which admit the above discussed situation, where 
the multidimensional Planck mass and the inverse size of the extra dimensions are both in the 
TeV energy range; the ADD scenario \21] in this case demands a very large number of extra 
dimensions. Moreover, a flaw of this approach is the assumption that the multidimensional 
background metric can be taken to be fiat, which means that the proper gravitational field 
of the brane can be neglected. The validity of the results obtained within this approach 
depends upon whether this approximation is good or not so good. In any case, studying 
the equations of motion for multidimensional gravity interacting with a brane of non-zero 
tension, it is not difficult to understand that if extra dimensions are compact, there should 
exist at least two branes, and the background metric must be essentially nonfiat. 

This situation is realized in the Randall-Sundrum model with two branes [22], - the RSI 
model, which is one of the most interesting brane world models. It is a consistent model 
based on an exact solution for gravity interacting with two branes in five-dimensional space- 
time. If our world is located on the negative tension brane, it is possible to explain the 
weakness of the gravitational interaction by the warp factor in the metric. A flaw of this 
model is the presence of a massless scalar mode, - the radion, which describes fluctuations of 
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the branes with respect to each other. As a consequence, one gets a scalar-tensor theory of 
gravity on the branes, the scalar component being described by the radion. It turns out that 
the coupling of the massless radion to matter on the negative tension brane contradicts the 
existing restrictions on the scalar component of the gravitational interaction, and in order for 
the model to be phenomenologically acceptable the radion must acquire a mass. The latter 
is equivalent to the stabilization of the brane separation distance, i.e. it must be defined by 
the model parameters. The models, where the interbrane distance is fixed in this way, are 
called stabilized models, unlike the unstabilized models, where the interbrane distance can 
be arbitrary. Below we will discuss the contact interactions of the Standard Model particles, 
which arise in the stabilized brane world model proposed in |23j . 

Let us denote the coordinates in five-dimensional space-time E = M4 x S^/Z2 by {x^^} = 
{x^, y}, M = 0, 1, 2, 3, 4, /i = 0, 1, 2, 3, the coordinate x* = y, —L < y < L parameterizing 
the fifth dimension. It forms the orbifold, which is realized as the circle of the circumference 
2L with the points y and —y identified. Correspondingly, the metric qmn and the scalar 
field (f) satisfy the orbifold symmetry conditions 

9f^uix, -y) = gf,uix, y), g^i{x, -y) = -g^A^x, y), (1) 
5-44(3;, -y) = 5(44(3;, y), 0(x, -y) = 0(x, y). 

The branes are located at the fixed points of the orbifold, y = and y = yb = L. 
The action of the stabilized brane world model can be written as 

S = j d^x j\yRV^+ J d'x Q/'^^a/ 09^0 - ^ (0) j (2) 



g\i{(f))d'^x + / ^/^{-\2{(^)) + LsM)d'^x. 

y=0 Jy=L 

Here ¥{(/)) is a bulk scalar field potential and Xi^2{4>) are quadratic brane scalar field po- 
tentials, g = detg^u, and g^u denotes the metric induced on the branes. The fields of the 
Standard Model are assumed to be located on the negative tension brane, Lsm denoting 
the Standard Model Lagrangian. A background solution in this theory, which preserves the 
Poincare invariance in any four- dimensional subspace y = const, looks like 

ds^ = e-^^'^y^^^dx^'dx" -dy^ = -iMN{y)dx^dx^, (3) 
<P{x,y) = (f){y), (4) 

the fields of the Standard Model being in the Higgs vacuum; it is also assumed that the 
vacuum value of the Higgs potential is equal to zero. For a special choice of the potentials 
[23] . the solutions for functions A{y),(p{y) are 

<P{y) = 0ie-"l^l, (5) 
A{y) = k{\y\-L) + ^{e-'-\y\-e-'-'^). 

The solution for A{y) is normalized so that the induced metric on the brane at ?/ = L is flat 
and the coordinates {x'^} are Galilean on this brane [TH [201 121] • The constants k,u, the 
boundary values of the scalar field 0i_2; the fundamental five-dimensional energy scale M, 
on a par with the coefficients of the quadratic brane potentials Ai^2(0); are the parameters of 
the model. When the former parameters are made dimensionless with the help of M, they 



4 



should be of the order 0(1), so that there is no hierarchical difference in the parameters. 
The separation distance in defined by the equation 



u 



(6) 



and, therefore, it is stabilized. 

As we have explained above, it is sufficient to treat the gravitational interaction per- 
turbatively to the linear order. To this end we represent the metric and the scalar field 
as 

gMN{x,y) = 'jMNiy) + -^j==hMN{x,y), (7) 

<P{x,y) = ^,e--\y\ + -J=f{x,y). (8) 

Substituting this representation into action ([2]) and keeping the terms up to the second 
order in l/\/2M^, we get the so called second variation Lagrangian which includes the 
Fierz-Pauli Lagrangian in the warped background, the terms describing the interaction of 
the linearized gravity in the background ([5]) with the branes, and an interaction Lagrangian 
determining the coupling to the fields of the Standard Model. In particular, the interaction 
Lagrangian is 

the energy-momentum tensor T'^" being canonically built from the Standard Model La- 
grangian LsM- 

In paper [2S] it was shown that the action built with the second variation Lagrangian 
can be diagonalized for any background solution and after the mode decomposition, which 
includes integration over dy, brought to the form 

^^ff = i E / {d'^b'^'^d^bl - + dx {d^^.d'^^u - fil^u^u) , (10) 

k=0 fc=l 

with mo = and the other masses of the four-dimensional tensor and scalar fields being 
defined by the background solutions A{ii) and 0(?/) and the model parameters. The free 
theory with this Lagrangian can be easily quantized, the propagators of the massive tensor 
and scalar particles are given by 



p pZ _ _|_ 



Uk ( \ 1 / Pt^Pp\ ( PuPa\ . 1 / Pt,Pa\ ( PvPp 



'^'^'"'^^^^ 2 V'''^^ ml ) ml 2 V'^'' ml J V'"' ml 

1 f PtiPv\ [ PpPa 



3 V'''^^ ml ) mi 
D\p) = , . (12) 
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The interaction of these four-dimensional fields with the fields of the Standard Model is 
defined by the interaction Lagrangian and looks like 



= --7^ ML)b%ix)T^'' + Y,^n{L)hl^{x)T^'' + \Y,9n{L)^n{x)Tli , (13) 

\ n=l n=l / 

ipn{y) and Qniy) being the wave functions of the modes in the extra dimension. Thus, 
the couplings are defined by the values of the wave functions on the brane. The latter are 
specified by the background solutions A{y) and (f){y) and the values of the model parameters. 
Since the field 6°,^(a;) describes the massless graviton, the coupling constant of this field to 
matter on our negative tension brane must coincide with the inverse Planck mass. The latter 
can be expressed in terms of the model parameters as [26] 



o_ f fk \ fk 



k2 



M?,^-e-,2r^|0(-,2.j-,^^.2.e-j|. . ^ ^, (14) 

7 denoting the incomplete gamma function and constant c = —kL — 6e^^"^. For u ^ 
(in fact, in this limit the background scalar field becomes constant, its fluctuations decouple 
from those of the gravitational field and the model goes to the unstabilized Randall- Sundrum 
model) this expression goes to the well known relation between the energy scales for an 
observer on the negative tension brane in the unstabilized RSI model 

2kL _ 1 

Ml, = M'^—, (15) 

which allows for a solution to the hierarchy problem of the gravitational interaction, if 
M ~ A; ~ ITeV and kL ~ 35 f2U[ 1^ . The problem of energy scales in the RSI model was 
discussed in detail in the paper [26]. In this paper the parameter space of the model was 
scanned and different scenarios for the fundamental energy scale and the KK excitations 
scale were studied. In particular, it was shown that it was possible to have the fundamental 
five-dimensional energy scale M of the order of 1-10 TeV with the masses of the tensor and 
the scalar KK excitations also being in the same energy range. The present day experimental 
data imply that this scenario is more likely, than the scenarios with the light radion. In this 
case the interactions of the Standard Model particles at the accessible energies due to the 
KK excitations of the tensor and the scalar fields can be very well approximated by a contact 
interaction, because we can drop the momentum dependence in the propagators. 

Integrating out the heavy tensor modes in the sum of Lagrangians (ITU]) . (IT^ induces the 
interaction of the Standard Model fields of the form 



vn>0 " 



= (16) 

111 

whereas integrating out the scalar modes induces the interaction of the form 
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The couphng constants eilF ^^^) '^^^ approximately esti- 

mated in the model as follows. For the stabilized Randall-Sundrum model presented in [25] , 
it was shown that it was more convenient to use parameters b = 0^/48M^, k = k — 2bu, and 
L instead of 0i^2 and k. It was also shown that for uL <ti 1 the metric of the stabilized model 
is similar to that of the unstabilized one with the inverse anti-de Sitter radius k instead of k, 
and it is possible to find analytical solutions for the wave functions of the tensor and scalar 
modes and their mass spectra. 

In particular, the spectrum of the tensor excitations is defined by Ji j = and in the 

approximation of an infinitely hard brane potential A2 we get ilJn\y=L = —\fk (see [SUES])- 
The sum over the tensor modes in (fT6|) can be estimated to be 



i)l(L) 2 1 0.1246 1.82 



1 ^ rn{L) 2 1 



2 ' 



8M3 ^ ml Ti^M^k ^ (1 + 4n)2 SM^ A^m? 

where we have introduced the coupling constant of the first KK resonance and its mass 
mi. 

mi = S.83k. 



It is worth mentioning that the contribution of the first KK resonance to this sum is exactly 
^^2^- One can also see that j- = — for relatively small n. 

In this parametrization, which is often used in the RSI model, the effective Lagrangian 
for this model takes the form 

Ujf = + = -^T^^A,.,,.T''^ (19) 

TT 1 

where 6 stands for the contribution of the scalar modes and will be calculated below. 

For the scalar sector the spectrum in this approximation is defined by [25j (where we 
have substituted mi instead of k) 

l + a + p.83 ■ _ 3.83 • (3.83- f.^^ 



mi J \ mi J mi \ mi 



with a = W^l + ^f^j + S^^^^T^ ~ 1-8 and for the wave functions we get 



^ '2 y/bu'^mi 



2 



We note that with given parameters mi, A^ and L, describing tensor sector of the model 
(one should take kL ^ 35 for the hierarchy problem to be solved), the spectrum of the scalar 
modes and their couplings to matter are also defined by the parameters u and b. Quite an 
interesting feature of the massive modes (both tensor and scalar) is that their masses and 
coupling constants in fact do not depend on the size of the extra dimension L, at least for 
relatively small n. 
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To estimate the corresponding sum over the scalar modes, we should specify the model 
parameters. Let us suppose that the lowest scalar mode, the radion, has the mass of the order 
of 2TeV. Such a situation can be realized if ~ 8TeV, mi ~ 3.83TeV^ (correspondingly, 
M ~ 2TeV, k ~ ITeV), u ~ O.OOSTel^, bu'^ ^ 0.28TeV'^. In this case the sum over the 
scalar modes in f|T8|) turns out to be 

« ^(0.341 + 0.002) «H11, 
^-^ ai mi mi 

n " 

where the first term corresponds to the contribution of the radion. Correspondingly, we find 
5 ^ 0.7. 

As we have mentioned above, this interaction Lagrangian leads to quite definite processes 
with the SM particles, which are determined by the structure of the energy-momentum 
tensor T'^'^. The latter is a sum of the energy- momentum tensors of the free SM fields 
and of contributions from the interaction terms, which are proportional to the SM coupling 
constants. The energy-momentum tensors of the free SM fields are quadratic in the fields 
and are explicitly given in Appendix A. 

One can easily see that for massless vector fields the trace of the energy-momentum tensor 
vanishes, and the scalar degrees of freedom do not contribute to the effective interaction. 
They can contribute to the effective interaction, if one takes into account the conformal 
anomaly of massless fields. The anomalous part of the energy-momentum tensor turns out 
to be 



^9) ( ^ d^d^ 



which gives the well-known expression for the anomalous trace of this tensor 

^~ 2g ' 

where b{g) is the beta function. The structure of this anomalous term in the energy- 
momentum tensor is such that the interaction due to the exchange of tensor particles fll6l) 
vanishes, and only the interaction due to the exchange of scalar particles (llSp remains. How- 
ever, this interaction is rather suppressed compared to the one due to the exchange of tensor 
particles, because the trace of the energy-momentum tensor is proportional to the particle 
mass that is much smaller than both mi and A^. And a possibility to observe the scalar 
component of the effective interaction may be due to the Higgs-radion mixing [27l [28] . 

Thus, in the lowest nonvanishing order in the SM coupling constants the effective La- 
grangian (fT9|) is a sum of four-particle effective operators (not only 4-fermions, but also 
2- fermions-2- vectors, 4- vector particles etc.). Experimental observation of production pro- 
cesses following from the effective Lagrangian (fT9!) or restrictions on their cross-sections 
allow one to estimate the multidimensional energy scale M, provided one gets a theoretical 
estimate for the product of the parameters mi and A^^ in fll9p . Their ratio may be estimated 
from the fact that the width of the first KK excitation must be smaller than its mass. 



3 Two body processes with KK gravitons 

As was mentioned, the lowest order effective Lagrangian in the SM couplings contains a 
sum of various four-particle (not only 4-fermions, but also 2-fermions-2-bosons, 4-bosons) 
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effective operators, which are gauge invariant with respect to the SM gauge group and lead 
to a well defined phenomenology. The Lagrangian involves only three free parameters A^, 
nil and S, where A^, mi parameterize the common overall coupling and 6 parameterizes 
the relative contribution of the scalar radion field (or fields as takes place in the stabilized 
RS model). In this paper we shall not present a detailed phenomenology, but rather point 
out some interesting aspects. In the leading order only the neutral currents of the same 
generation SM fields are involved. These new interactions do not lead to additional decay 
modes. Possible new decays of the SM particles from the effective Lagrangian may only be 
present in the next order in the SM couplings, when charged currents appear in the SM 
energy-momentum tensor. Also new effective 4-particle operators following from the SM 
energy-momentum tensor obviously do not lead to flavor changing neutral currents. In the 
tree level approximation there are several processes following from the effective Lagrangian, 
which appear only at loop level in the SM such as gg I'^l' , gg ZZ(W^W~), e+e" — >• gg, 
11 ~^ 99 6tc. In Appendix B analytical expressions for the total and differential cross 
sections for the processes gg l^l^, gg ZZ{W^W~), qq l^l'^, qq ZZ{W~^W~), 
e~^e~ f f, e~^e~ gg, 77 f f, 77 — >• gg are presented. For completeness we keep 
nonzero masses of the final state particles. In the case of massless fermions formulas for the 
total and differential cross sections for the Drell-Yan processes gg — > and qq 
are in complete agreement with [161 ISH EQ]. Formulas (!38l) - (HTl) for processes gg ZZ and 
qq ZZ that take into account scalar KK modes, massive final states and the interference 
with the SM amplitudes are presented here for the first time. In the cases where colliding 
gluons produce massive final particles, there is also a scalar radion contribution, which is 
proportional to the parameter 5^ of the order of 1 and to the trace anomaly coefficient 
{h{gs) /2gsY . We give this contribution in formulas, although numerically it is about 100 
times smaller than the corresponding tensor contribution. 

Below we will perform numerical simulations for the Drell-Yan process because this chan- 
nel is most sensitive for new physics. Detailed simulations for other channels will be made 
in a further study. 

Symbolic and numerical computations, including simulations of the SM background in a 
thought experiment for Tevatron and LHC, have been performed by means of the version 
of the CompHEP [31] package realized on the basis of the FORM [32] symbolic program. 
The Feynman rules following from the effective Lagrangian have been implemented into this 
version of the CompHEP. Such an implementation allows one to use the code for event 
generation and to perform analysis in future more realistic studies. 

Qualitatively, the situation from the phenomenological point of view is similar to that 
appearing in the ADD scenario and worked out by J. Hewett in [16]. The correspondence 
between the parameters used in our study and in [16] is the following: 

A>f _ IMl 
0.91 ~ 4~' 

As shown for the RSI model in [17] the exchange of a tower of the KK gravitons in the energy 
range below the KK production threshold, similar to the ADD case, leads to an increase of 
the invariant mass tail of produced particles. For the Drell-Yan process it is demonstrated in 
Figs, dl [21 The process gg contributes to the Drell-Yan process and it was included 

in our numerical simulations. As shown on Fig [51 this contribution is very significant for the 
LHC As will be demonstrated below, even in the case when the first KK resonance lies in 
the energy range accessible for a detection one should take into account the contribution 
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from all the other KK states. One should stress that in the ADD scenario, in addition to 
deviations from the SM prediction for the processes such as lepton pair production, there 
should also exist processes with the KK tower radiation off. The latter processes do not take 
place in the RS model. 

Using the standard analysis and taking into account the expectations for system- 
atic uncertainties (detector smearing, electroweak, QCD scale, Parton Distribution Func- 
tion(PDF)) and statistical uncertainties of the SM dilepton invariant mass shape (see exper- 
imental data [33] for the Tevatron and Monte Carlo simulations [34j for the LHC), we obtain 
the current Tevatron limit for the coupling parameter at 95% CL and estimate expected 
experimental limits for this parameter (Table [1]) that may be reached at the Tevatron for 
higher luminosities and for various luminosities at the LHC. We used CompHEP for the 
calculations of the SM center values in thought experiment. 



Table 1: Experimental limits for the coupling parameter at 95% CL that may be reached at 
the Tevatron and the LHC using Drell-Yan process for some values of integrated luminosity 





TEVATRON {^s = 1.96 TeV) 



LHC (v^ = 14 TeV) 



L, fb-' 


at 95% CL, TeV-^ 




at 95% CL, TeV-^ 


1 


1.185 


10 


0.238 -10-2 


2 


0.995 


20 


0.203 -10^2 


3 


0.900 


30 


0.184 -10-2 


5 


0.790 


50 


0.164 -10-2 


10 


0.664 


100 


0.140 -10^2 



The Tevarton limit for 1/6 ^ of integrated luminosity expressed in terms of parameter 
^GRW introduced in [35j 



m: 



GRW 



0.91 



2ti A-lml^ 

gives M^^^(l/6^^) = 1.52 TeV, which is in a good agreement with the corresponding limit 
from the cited experimental paper |33j . 

The last string of Table [T] contains limits corresponding to the highest value of collider 
luminosity: 



Tevatron{10fb 



0.91 

A^mf 



xTeV^KOM, LHC{100fb 



-1^ 



0.91 

A^mf 



xTeV^ < 0.0014. (21) 



Figures [3] and m demonstrate distributions corresponding to values (1^ . These limits may 
be used for estimating the lowest value of parameter A^^ from a requirement that the width 
of a resonance be smaller than its mass: Fi < rrii/^, where ^ is some number, ^ > 1. Using 
limits fl2T]) and the equation for the total graviton width (|5T|) ^™^^ |^ < we get 



Tevatron : A^ > 0.61 ■ ^^/^ TeV, LHC : A^ > 2.82 ■ ^^/^ TeV, ^ > 1. 



(22) 



One of the effects in searches for KK resonances below the production threshold of the 
first state is an enhancement of the effective coupling due to KK summation in comparison 
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Figure 1: Dilepton invariant mass distri- 
bution for parameter x TeK^=0.66 

(dashed-dotted line), 1.82 (dashed line), 4 
(dotted line) for the Tevatron 



Figure 2: Dilepton invariant mass distri- 
bution for parameter x TeV^=0.0014 

(dashed-dotted line), 0.0046 (dashed line), 
0.01 (dotted line) for the LHC 




Figure 3: Dilepton invariant mass distribu- Figure 4: Dilepton invariant mass dis- 
tion for 95% CL parameter ^^^^ xTe\^^=0. 66 tribution for 95% CL parameter x 

for the Tevatron (L = 10/6"^)' Te1/^=0.0014 for the LHC (L = 100/6"^) 
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Figure 5: Dilepton invariant mass distribu- 
tion for 95% CL coupling parameter for the 
LHC. The dashed-dotted hne corresponds to 
SM + KK sum {gg — > fi'^fi'), the dashed hne 
corresponds to SM + KK sum {pp —>■ ^'^ 
process. 




Figure 6: Dilepton invariant mass distribu- 
tion from the SM (solid line) and from the 
SM plus sum of KK modes including the first 
KK resonance with M^es = 3.83 TeV, Tj-es = 
0.08 TeV, A^ = 8 TeV (dashed line) for the 
LHC 




Figure 7: The normalized dilepton invariant Figure 8: The normalized dilepton invari- 
mass distribution from the first KK resonance ant mass distribution from the sum of KK 
plus the sum of KK tower states starting from tower states starting from the first KK mode 
the second mode (solid line) and from the first (solid line) and from the first KK mode only 
KK resonance only (dashed line) for M^es = (dashed line) for M^es = 10 TeV, F^es = 
3.83 TeV, Tres = 0.08 TeV, K = 8 TeV for 0.5 TeV, = 14 TeV for the LHC 
the LHC 
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to the first mode contribution below the threshold only. For the considered case of the 
stabilized RS model one has a factor 



This leads to an increase by 1.82^ ^ 3.3 times in the production rate (for the case of one fiat 
extra dimension this factor is 1.64 => X] ^ = being numerically close to the warped 
case) . 

To illustrate changes in distributions due to KK tower contributions we run simulations 
for two parameter points with the first KK resonance being in and out of directly detectable 
regions. The first point (mi = 3.83 TeV, A^r = 8 TeV, Fi = 0.08 TeV) was already discussed 
in Section 3. Such an RS resonance (see Fig. [6]) is close to the direct reach limits expected 
for the LHC [17J. For the second point (mi = 10 TeV, = 14 TeV, Fi = 0.5 TeV) the 
mass of the first KK excitation is close to the collider energy limit, and it is not directly 
observable. For both points we can use the low-energy effective Lagrangian approach. The 
effective Lagrangian allows us in both cases to sum up the contributions from all the KK 
modes or from all except the first one, and in this way to take into account their influence 
on the background tail. As one can see from Fig. [7] and Fig. [H], the additional substrate 
from the KK tower increases the production rate more than 3 times in the invariant mass 
region below the resonance mass. The situation is significantly different above the resonance, 
where in addition to the resonance pike there is an area with a minimum due to a destructive 
interference between the first KK resonance and the remaining KK tower contribution. This 
local minimum takes place at the value of invariant mass M^m ~ 1.5mi. The growth of the 
invariant mass after the minimum is strongly suppressed by parton distribution functions 
leading to an additional bump in the invariant mass shape. But this bump is unlikely to be 
visible in the experiment on top of the SM background as shown in Figj6l 

In conclusion of this part, one should stress that in order to perform correct searches for 
KK resonances not only interferences with the SM if nonvanishing and computed NLO QCD 
corrections should be included into corresponding generators, but also the influence of 
those KK states, which are not reachable directly. 

4 Conclusion 

In the present paper we have continued studies of the effects that appear in theories with 
extra space dimensions due to the exchange of directly unaccessible KK modes of the fields 
propagating in the bulk when the fundamental energy scale that defines these masses is larger 
than the typical collision energies 1 — 10 TeV. We have derived the effective Lagrangian re- 
sulting from the exchange of these KK modes in the stabilized RS model. Such a Lagrangian 
has a simple structure of a product of two currents, corresponding to the zero mode of a bulk 
field, multiplied by an effective coupling constant of a dimension depending on the spin of 
the bulk field. The exact value of this coupling is model dependent and has to be calculated 
for each model separately. In fact, this structure is reminiscent of the Fermi interaction. 
Obviously, if the bulk field is the gravity field, then the corresponding current is just the 
standard energy-momentum tensor of the SM fields (see Appendix A). A delicate nuance in 
the case of the stabilized RSI model is that the multidimensional metric has both tensor and 
scalar massive degrees of freedom, the latter coupled to the trace of the energy-momentum 
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tensor or to its anomaly in the case of massless SM fields. It is shown explicitly that the 
contributions of the scalar modes are much smaller than those of the tensor modes. 

The Feynman rules for the new effective 4-particle interaction vertices of the Standard 
model fields have been incorporated into the CompHEP computer program, and explicit 
formulas for the differential and total cross sections for a number of processes, generated by 
these effective Lagrangians, are presented in Appendix B. For completeness we have included 
into the symbolic formulas the final state particle masses and contributions of virtual graviton 
and radion KK modes. 

Using these results, we have calculated and plotted the contributions due to the multidi- 
mensional gravity to the Drell-Yan production processes for the Tevatron and LHC energies 
for a number of the stabilized RS model parameter points. It was clearly demonstrated that 
an enhancement of the effective coupling due to KK summation in comparison to the first 
mode contribution only leads to an increase of the collider potential to probe the first mode 
mass below the threshold of its production, the mass being significantly larger than the 
collision energies. In the case, where the first mode mass is in the accessible energy range, 
but all the other modes are out of this range, summation of all the modes contributions 
starting from the second one leads to a significant change of the shape of the Breit-Wigner 
distribution. The latter case occurs for the considered explicit example in the stabilized 
Randall-Sundrum model. 

We have calculated the effective couplings for this case. To this end, we have used both 
analytical results and numerical estimates for the wave functions of tensor and scalar modes 
and performed an approximate summation of the series of inverse mass squared of the KK 
excitations. For a choice of the parameters of the model, where the effective energy scale 
M ~ 1 TeV and the other parameters are such that the metric of the stabilized RS model 
with the radion mass of the order of 2 TeV is similar to that of the unstabilized one, we found 
explicit formula f|T9l) for the effective Lagrangian. For this particular choice of the parameters 
the first KK graviton and the radion masses are beyond the energy range directly accessible 
at the Tevatron and within this range at the LHC. However the value of the effective coupling 
in (fT9|) is too small for deviations from the SM tails to be observed and to probe the masses 
below the production threshold at the Tevatron energies. For the LHC the first KK and 
the radion masses are within the directly accessible region. In this case in order to perform 
correct searches for the tensor resonances and to model the distribution tails one should sum 
up the contributions from all the other KK modes and take into account their interference 
with the resonances. It is worth pointing out that this summation of the contributions of the 
massive tensor modes is also needed for modeling the background in searches for the light 
radion, whose coupling to the SM fields turns out to be larger, than in the case of the heavy 
radion. 
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5 Appendix A: Energy- momentum tensors for free SM 
fields 

The Lagrangian and the energy-momentum tensor for fermions: 

= 1 - 9/,%^*) - m*** (23) 



= \ (^7m^<^* + *7.5/.* - 5.*7/.* - 5m*7.*) - (24) 

T*^ = (*7''5;.* - 5^*7''*) + 4m*** (25) 

The Lagrangian and the energy- momentum tensor for massive vector bosons (Z-boson): 

Lz^--Z,,Z>^^ + ^Z^Z^ (26) 

T^, = -Z,,Z,,g^^ + mlZ.Z, + r^,, Q^pa^"'^ - ^^'^p) (27) 

T'l = -m\Z^Z, (28) 
The Lagrangian and the energy-momentum tensor for complex vector bosons ( W-bosons) : 

W = -^^t^"^-^^ + m%,W;W-^ (29) 
= -W+W-^g'^'^ - W:^;W;,g^'^ + w?^ {W^W; + W^+W^") + (30) 

T^l = -2ml,W;W-' (31) 

The Lagrangian and the energy-momentum tensor for massless vector bosons (the photon 
and the gluons): 

La = -\f,.F^'' (32) 

T^t = -F.pF^aQ'"' + V,u\Fp.FP'' (33) 
= (34) 
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The Lagrangian and the energy-momentum tensor for the scalar field (the Higgs field in 
the unitary gauge): 

= Id^'^d,^ - (35) 

^^') (36) 

(37) 

6 Appendix B: Partonic total and differential cross sec- 
tions for 2^2 processes. 

Processes pp Z^Z^: 



da, 



gg^zz 
dz 



167r • 16 



(3 



36' 



%) 
2g 



(38) 



K 



^gg-^ZZ — 



167r-8 



p 



(39) 



da, 



qq-^ZZ 

dz 



IGtt ■ 96 
16-48 

TT • a^C2 

16-48 



(3 [A- 2(3\z'' + 1) + 3/3V(l - z^)] 



(40) 



+ /3^(2/3;. + /3^-3) + e2'/? 



. 6 

- ^ (Ci - ^)' ■ e 



+ 



^qq-^ZZ 



K 



167r • 12 

KCK • Ci 

16-48 

TT - a;^C2 
16-24 



2,(1-/3^ 



(5 - /3^) 



In 



1 + /3 



1-/5 



s+ 



1 + + 4 



1-/3^ 
I + /32 



In 



1+P 



1-/9 



2/3 



(41) 
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Processes pp l^l : 



dz ~ 167r • 64 



[2{i - z^) -13'' {I -z'^f] 



K 



(^gg^l+l- 



16n ■ 24 



5 



da, 



qq-*l+l- 

dz 



K 



(^qq^l+l- 



-f3^ [iz^ + l) + Aph^{z^-l)] s'+ 



IGtt ■32-3 

^ [zP^{l -P^ + z^(3^) {2Q, - 2viVgXi) + (1 - 3z'')(3'aia,xi] s+ 



2-3 



[(2/3 - + fi^z") {Ql - 2xiQqViVg + X2{a^ + ^;f)(aj + ^J)) + 



2p^z(-2xiQq + Ax2ViVg)aiag + 2x2(af + vf)al(P^ - 



1671 ■ 36 
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(2/3 + -P') {Ql - 2xiQqViv, + X2(af + vf)(al + v^^)) + 



2x2{af + vf)al{P' - P)] ^ 



Processes for linear colliders: 



da, 



e+e~— >g+q'~ 

dz 



= 9- 



da. 



qq—*e~>~e~ 



dz 



dae+e--^gg , 4\ ^3 
; — — (1 — Z )S 

dz 167r • 2 ^ ' 



^e+e -^gg 



K 



Att ■ 5 



d^2p9l = -^{1 + 6^2 + /)s^3 
dz lOTT ■ 2 
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TT • O 

Total width for the KK graviton resonance 



Ml 



AE ■ 47r ■ 5 



9 1 /13 14M| M|\ ^ 2 /13 MM^ , ^ 



21 /3_lilg_ M^- 
8 V 8 2 M2 Mi ' ^* 



n n 



JW|_97 
A2 • 47r 80 



where: 

S — X1X2S, 



AM} 

Pf — \j 1 r^, Mf standing for the mass of the final state particle, 

z = cos 9, 9 denoting the partonic center of mass scattering angle, 
9w being the Weinberg angle, 
0.91 

K — — — 7T, with A^ being a coupling constant of the first KK resonance and its mass mi, 
A^mf 

ai = -1, 

vi — —1 + 4sin^ 9w, 

dq — '2T3q, 

Vq = 273^ - AQgSin'^9w, 

Qq denoting the charge of the quark, 

T^q = ±- for up/ down — type quarks, 
2/5 ' 2/3 ' 



a\ + vl 6a^v^ -^a^^vl 



Xi 

X2 



sin^ cos2 9w ' sin"^ cos^ 9w 

1 s(s-M|) 



16 sin^ 9w cos2 - + M|r| ' 

1 s2 



256 sin^ 9w cos^ (s - M|)2 + M|r| ' 
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